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Abstract. 

The aim of this paper is to revisit the calculation of atom-surface quantum friction 
in the quantum field theory formulation put forward by Barton [Neiv }. Phys. 12 
(2010) 113045]. We show that the power dissipated into field excitations and the 
associated friction force depend on how the atom is boosted from being initially at 
rest to a configuration in which it is moving at constant velocity ( v ) parallel to the 
planar interface. In addition, we point out that there is a subtle cancellation between 
the one-photon and part of the two-photon dissipating power, resulting in a leading 
order contribution to the frictional power which goes as v A . These results are also 
confirmed by an alternative calculation of the average radiation force, which scales as 
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1. Introduction 

The interaction of moving objects with light has been in the focus of physics even 
before Einstein's annus mirabilis and his fundamental papers about special relativity. 
A seminal contribution in this context is Einstein's derivation of Planck's blackbody 
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radiation law (HO: it has brought upon us not only the concepts of spontaneous and 
stimulated emission. Also the momentum exchange between atoms and photons, 
and the corresponding friction and diffusion have been shown to provide the 
physical picture for the thermalization of the velocity distribution of an atomic gas, 
decades before the advent of laser cooling techniques |3j|. Radiative friction (without 
lasers) is the process where a moving atom comes to rest in the preferred frame set 
by the blackbody radiation field |4). (Motion relative to the frame of the cosmic 
microwave background, for example, can indeed be detected by the anisotropy in the 
apparent temperature 0 .) Quantum friction is the theorists' variant of this problem, 
when the temperature is set to zero. Velocity-dependent (or drag) forces only appear 
when true relative motion is defined by the presence of another object. In this paper, 
we consider the simple case of an atom (or molecule) near a macroscopic half-space 
filled with metallic material. The distance between the atom and the metal surface is 
also macroscopic (at least a few nm) in the sense that electronic overlap is negligible. 
In this regime of distances, it is valid to use a local approximation for the optical 
response of the surface, its permittivity depending only on frequency. 

It is instructive to draft a short summary of the long series of works dealing with 
the problem of quantum friction on an atom moving at constant velocity parallel to 
the vacuum-metal interface. We will restrict ourselves to works that mainly used the 
local approximation for the optical response, i.e. those that considered macroscopic 
distances in the sense defined above. It is interesting to note that various authors 
obtained quite different results for this drag force, differing both in their dependence 
on velocity and with atom-surface separation. Unfortunately, most of these works do 
not critically discuss the others nor attempt to clarify the origins of the differences. 
One of the earliest works on the problem was undertaken by Mahanty flU, who 
computed the velocity dependence of the drag force on a moving molecule. It was 
found that the quantum friction force scales as vz~ 5 for small velocity v and large 
separation z between molecule and surface. However, this calculation was criticized 
by various authors since it predicts a non-zero quantum friction even for a perfectly 
reflecting surface (which lacks the possibility of referencing relative motion, indeed). 
Another series of papers also obtained a linear dependence of quantum friction on 
velocity. Schaich and Harris |7J computed dynamic corrections to van der Waals 
potentials for a neutral molecule moving above a metallic plate, and modeled the 
molecule as a dipole oscillating normal to the surface. The resulting friction force is 
again linear in velocity but with a different asymptotic large distance dependence as 
z l0 . More recently, Scheel and Buhmann 0 have considered a multi-level atom 
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moving at constant velocity, and employed a master equation approach to solve 
for the atom dynamics in the Markov approximation. They again found a linear 
dependence on velocity and a z~ 8 scaling in the near-field. These same scalings (with 
slightly different numerical pre-factors) were obtained by Barton 0 in a harmonic 
oscillator model for the atom, where the friction force is computed in time-dependent 
perturbation theory from the power dissipated into pairs of plasmons. Hoye and 
Brevik have put forward an approach to quantum friction very similar to Barton's, 
and used it to compute the friction force between two atoms IfTOl or two plates IfTTH 
and compared it to Barton's results for these particular systems. 

In contrast to all the above works in the literature, various other authors have 
obtained a vanishing contribution to the atom-surface friction force linear in velocity. 
For example, Tomassone and Widom Ifl2l computed the finite temperature friction 
force on molecules moving near metals using the image charge approach, and 
obtained a vanishing linear-in-v quantum friction in the limit of zero temperature. 
The same conclusion was reached by Volokitin and Persson Itl3ll . who employed 
fluctuation electrodynamics to compute the Lorentz force on a moving dipole, by 
Dedkov and Kyasov lTl4ll . who used the equilibrium fluctuation-dissipation theorem 
to evaluate the dipole and field correlation functions, and by Golyk, Kruger and 
Kardar fllSfl. who evaluated the force using linear response relations in fluctuation 
electrodynamics. Another series of papers confirmed these results, and derived 
the first non-vanishing contribution to the quantum friction force that scales as v 3 . 
These include the works of Dedkov and Kyasov Ifl6l . who extended their previous 
calculations to capture the nonlinear dependence of quantum friction on velocity, 
Pieplow and Henkel ||T7|. who used equilibrium fluctuation electrodynamics to 
derive a relativistically covariant formulation for the friction force, and Intravaia, 
Behunin and Dalvit ffl8l , who calculated the atom-surface drag force by generalizing 
fluctuation-dissipation relations to the non-equilibrium stationary state defined by a 
constant velocity. 

One of the goals of this paper is to revisit the calculation of quantum friction in 
probably one of the simplest and cleanest formulations of the problem put forward 
by Barton 0. Within this approach, the zero-temperature friction force is computed 
by solving the joint atom+field/matter dynamics in time-dependent perturbation 
theory, starting from an initial state in which the atom and the field/matter 
subsystems are both in their ('bare') ground states. As emphasized by Barton, this 
perturbation theory has no need to rely on assumptions related to correlation times, 
linear response, or local thermodynamic equilibrium which are implicit in many 
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calculations performed with the toolbox of fluctuation electrodynamics. One also 
does not require fluctuation-dissipation relations. The challenge of this approach 
is that the dissipation in the atomic system is purely radiative and is generated 
self-consistently in the perturbation series. This is in sharp contrast to a field 
theory like the one reported by Volokitin and Persson Ifl9l where the basic two- 
point functions for atomic variables are constructed by a re-summation procedure 
including radiative damping. For simplicity, we will restrict ourselves to the near¬ 
field regime, where quantum friction is expected to be enhanced. We demonstrate 
in particular that the power dissipated into field excitations and the associated 
friction force depends on how the atom is boosted from being initially at rest to 
a configuration in which it is moving at constant velocity parallel to the planar 
interface. 

The paper is organized as follows. Secj2] reviews the building blocks of the 
quantum field theory for the atom-field interaction and gives the time-dependent 
state including amplitudes for one- and two-photon processes. In Secj3]we use these 
results to calculate the frictional power and force in the case where the velocity of the 
particle is constant for all times. Although this obviously requires an external energy 
supply to compensate for the frictional loss, the description is actually simpler, and 
one recovers some of the results presented in Ref.|[9|. It is shown in particular that 
the 0(v 4 ) contribution to the power of two-photon emission found in Ref.f9l (called 
there Pa) can be explained in terms of this special trajectory. Sec j3] also provides an 
alternative picture where the expectation value of the force operator is computed in 
the time-dependent state. Its stationary value at long times is found to scale with 
the velocity like ~ v 3 . Sec jl] contains the main results of this paper. The calculation 
of the radiated power is generalized to more realistic trajectories where the atom 
starts at rest and is accelerated to a constant final velocity. We discuss the role 
of the finite duration of the acceleration and show that: (i) the results presented 
in Ref.H depend of the specific choice of the atom's trajectory; (ii) that the power 
bookkeeping in Ref.[9]] is incomplete and needs to be complemented with the power 
needed to create the excited state. If this is done, we again find a frictional force that 
scales as ~ v 3 with the velocity. Sec j5] provides a review of two approaches |8lIT8| 
that describe quantum friction within the framework of fluctuation electrodynamics. 
Some technical material is relegated to the appendices. 
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Our discussion is based on Refs.|9ll20| where one of the simplest field theories for 
atom-photon interactions is developed. The physical situation is sketched in Fig|l] 
(left): a point-like atom moves at constant velocity v parallel to a half-space that 
responds linearly to the electromagnetic field. The distance z — za of the atom, kept 
fixed, is taken much smaller than the relevant wavelengths (non-retarded regime) so 
that the field can be described by an electric potential <3> [Eq.Q below]. The half¬ 
space is absorbing light, broadening the surface plasmon resonance. For simplicity, 
we still call 'photons' the elementary excitations of the field, although 'plasmon- 
polariton' or 'medium-assisted polariton' would be more appropriate names. The 
atom is described by a few low-lying states (FigJlJ right), and its position follows a 
prescribed trajectory r(£). Our goal is to calculate the radiative friction force F and 
the frictional power P = —r ■ F that must be supplied by the external agency that 
keeps the atom on its path. 



Figure 1 . (left) Sketch of the considered geometry, (right) Atomic energy levels. 


2.1. Relevant states and observables 


The lowest quantum states of the atom are taken by analogy to the Is and 2p level 
of the hydrogen atom: they are denoted | g) for the Is state, the three degenerate 
2pj states are written \fj). The unit vector fj is taken from a set { fj} forming an 
orthonormal basis that we may assume real without loss of generality. The Bohr 
transition frequency between the levels is Ft. We focus in this paper on transitions 
among these energy levels only and mention briefly where additional states would 
appear. The nonzero matrix elements of the electric dipole operator D(t) in the 
interaction picture are 


(9\D(t)\v) =vde 


—i 


(v\D{t)\g) = fjde m , 


(1) 
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where the transition dipole matrix element d is the basic coupling constant of the field 
theory. It determines, for example, the static polarizability a = (2d 2 )/(HQ) [Eq.(2.2) 
of Ref.U9j]. Using rotational symmetry, one also has the identity Y.fj r h r h — $ij when 
the three excited states are summed over. 

The atom+field coupling (summation over double indices is assumed hereafter) 

V(t) = -Di{t)Ei(r(t),t) = Did&{r(t),t) (2) 

is explicitly time-dependent via the atomic trajectory r(t). We will also often use 
the notation r (t) for the path in the xy-plane parallel to the surface placed at z = 0. 
In our approach the force acting on the atom parallel to this plane is given by the 
operator ETil 

F(t) = D l (t)VE l (r(t),t). (3) 


More general situations would include higher-order multipole moments of the 
atomic charge and current distribution (magnetic dipole, electric quadrupole ...) and 
the time derivative of the electromagnetic momentum d x B. The latter includes 
the so-called Rontgen interaction that takes into account the transformation of the 
electromagnetic field into the frame co-moving with the atom. This interaction is 
relevant at larger (retarded) distances |j8|. 

The field operator is expanded in a plane-wave basis of elementary excitations 
(photons) and evolves freely according to 


OO 

<£(r, t) = Jd 2 k j dcu (d ka , d kw exp(ik • r - i cot) + h.c.) . (4) 

o 

The bosonic operators satisfy the commutation relation [d kaJ , d k ,^,] = <5(k — k')5(a; —a/) 
and the 'one-photon amplitudes' are given by [Eq.(2.4) of 1191 1 


— 


yJujTuj 2 /2 rY 


—kz 


uj 2 + icejT — (jjl V 27 r 2 k 


(5) 


where the frequencies u p , c us, and T parametrize the dielectric function e(u) of the 
half-space. We note in particular the relation 


I0kj 2 = 


h e 2kz _ s U e 2kz , e(u) — 1 


Irn R(uj) = 


Im 


( 6 ) 


2tt 2 k y ’ 2i r 2 k e(o;) + 1 ’ 

where R(u) is the non-retarded reflection coefficient of the surface. The frequency ojs 
gives the surface plasmon resonance and T its broadening (complex pole of R(u)). 
With this expansion for the photon field, the force operator (j3j, for example, takes the 
form 


F(f) = / d 2 k / dta ( k (D(t) ■ k)f exp(ik ■ r(t) — i ut) + h.c. 


(7) 


o 
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where the three-dimensional wave vector k = (k, i k) is in fact complex with k ■ k = 0 
and k ■ k* = 2k 2 . Note that this is a 'skew' operator that connects quantum states with 
different photon numbers (for the field) and different energy levels (for the atom). 

Since part of our focus will be on the power radiated into photons and pairs of 
photons, let us introduce 


Pi = lim Y' /d 3 nh(Fl + 


UJ 




P 2 = ^ Hm Jd 3 Ki J d 3 K 2 Ti(uji + uj 2 ) 


\(g,KiK 2 \^(t))\ r ^ 


( 8 ) 

(9) 


where \{ff, kco|(t))| 2 and \{g, kitaik 2 a; 2 |'t(f))| 2 are the probabilities of finding the 
atom in an excited state and one and and the two-photon, respectively (see Secs |2.2| 
and [4]). We also used the compact label n = ku and the factor \ accounts for double 
counting the symmetric two-photon states. The long-time limit is to be understood 
within time-dependent perturbation theory: t is typically not longer than a fraction 
of the relevant life times. The two-photon power has been calculated in Ref.|9l; we 
review the evaluation of the integrals in Appendix A[ This calculation is generalized 
in Sec|4j where also a partial cancellation between Pi and P 2 is found. 


2.2. Atom+field states 

In our perturbative approach the initial state is given by the tensor product of atomic 
ground state and zero photons, \g, vac), while the interaction is represented via the 
operator V(t). An expansion up to the third order in the coupling constant d of the 
atom+field state |\l/(£)) leads to 

W)> = (l + Co 2) (t))|s,vac) + Ysf d 3 K(c ( i\t) +c ( i\t)MY 

V 

+ \J' d 3fi + d 3/t 2 4 2) (0 \g,KiK 2 ) + ... (10) 

where cfft) denotes the transition amplitudes for states with n photons in the 
p th perturbative order and can be obtained by using the standard techniques of 
perturbation theory. The relevant matrix elements are given by 

{g,yac\V(t)\fj, k) — id (ff ■ fc)0 K exp[—i(fl + u)t + ik • r(f)] , (11) 

(f},K\V(t)\g, K 1 K 2 ) = i d(rf- kf)f Kl e l{n ~ Ul)t+lkl ' r ^S(K - n 2 ) 

+ i d(ff- kf)f K2 e i(«-«2)t+ik2T(t) ( j( /c _ Kl ) ; (12) 

{if, vac |V"(f) | g, n) — id {if- fc)0 re exp[i(Q — ui)t + ik • r(f)] . (13) 







Friction forces on atoms after acceleration 


8 



g vac 


g kik 2 


Figure 2. (Color online) Schematic overview of the relevant states for the atom+field 
theory. The atomic states are labelled g and e (actually, e is threefold degenerate). The 
notation n = ku collects the quantum numbers of the field states, vac is the vacuum 
state. Two-photon states are denoted k+k 2 . To the left of the vertical dashed line are 
states that are not directly coupled to the ground state by the atom+field interaction, 
illustrated by the dotted blue lines. The thick black lines indicate the position of the 
surface plasmon resonance wg, the shaded grey area illustrates its broadening (not to 
scale) over a range T (intrinsic) and due to the Doppler shift k • v. The Doppler shift 
is responsible for the appearance of negative frequencies, intuitively understood as 
measured in the frame co-moving with the atom. The thin dash-dotted line illustrates 
the two-photon process gv ac ++ gn\ k 2 - 


It is important to note that the matrix elements depend on the detail of the path r(f). 
Let consider first the simple case of a constant velocity, i.e. r (t) = v/. (Corrections 
arising from a realistic trajectory including an acceleration stage ('launch') are 
discussed in Secj4j) In this case time-dependent perturbation theory leads to 

i dfq-kyfi 


r {1) 

c i 


c (2) 

c 2 


(*) = 

(0 = 


, . . v exp[i(n + u/)f] 

_ d 2 (k 1 ■ k 2 )*(j)* K1 (j)* K2 exp[i(u;( + u' 2 )t] 

Ti + co 2 — iA) 

1 +{10 2 } 


(14) 


x 


.fl + cuj - iA ' ^ J (15) 

where the positive infinitesimal A ensures that the atom-field interaction is 
adiabatically switched on in the past. The scalar product (k 1 ■ k 2 )* = k x • k 2 — kik 2 
arises from the summation over the three excited states [if) and we have used the 
shorthand u' = oj — k • v for the Doppler-shifted frequency (as 'seen' by the moving 
atom). 

The second-order correction Cq ( t ) needs a special handling because it involves 
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the energy shift of the state \g, vac) and the rate for the process | g, vac) —> \if, k): 


- ; 




K 


d 2 \f) ■ k\ 2 \f K 


fr(fl + uj — k • v — iA) 


dr e 


/ _Ar 


(16) 


The second integral integral formally evaluates to A 'e A/ « t + A 1 for A —>■ 0 so that 
we identify the frequency shift and inverse lifetime from 1 + Cq 2) (t) ~ exp(—i 8Egt/Ti — 
r ) g t/2). This yields 


-i SEg 


hfh = • V f j 3 d 2 |rf- k\-\<f) K \ 2 

2 + to — k ■ v — iA) 


(17) 


The level shift 5E g arising from Eq.([l7j) has been worked out by Barton, Eq.(3.2) of 
Ref.|9]|. For the lifetime, we get the 'Golden Rule' result ||30| 

Ott _ r ^ 

7 j = ^E J d 3 n\(fj, K\V\g,v&c)\ 2 5(Fl + a; - k ■ v) , (18) 


where the atomic motion leads to the Doppler shift of the final photon frequency 
This integral is, however, exponentially small for reasonable parameters, as can be 
seen as follows. The sum over the excited states gives J2if\v • k\ 2 = k ■ k* = 2k 2 . 
Inserting the coupling strength \ f K \ - from Eq.Q, we get 

/V/O /* 

7 „ = — / d 3 K ke~ 2kz Im R(oj) 5(Q + u — k • v) . (19) 

The resonance condition 0 = fl + ca — k v jointly with uj > 0 limits the domain for the 
^-integration to k • v > Q. We then have k > Fl/v, and the exponential gives a scaling 
of this integral proportional to - lz ! v _ We follow here the same strategy as Ref.|9l 
and neglect contributions that show such an exponential scaling, assuming that the 
velocity is small enough: v Flz. [For a lithium beam at 10 keV and 10 nm distance, 
v/Vtz ~ 0.02.] The physical interpretation of this process is the following [see also 
Ref. [[311 in this issue]: the atomic motion leads to an anomalous Doppler shift (a/ < 0 
in the co-moving frame, while uj > 0) that makes the 'spontaneous excitation' of 
the ground state possible, similar to Cherenkov radiation H32l[33[. The rate for this 
process is, however, extremely slow because of the threshold set by the atomic Bohr 
frequency, uj' = —Q. 

For the third-order correction to the one-photon process one gets: 

t 

Ci\t) = ~ ff J d3fi: i d3/t 2 J (ff, n\V(t')\g, 

— OO 

t 

-jr J dt' (ff,n\V(t')\g,vac)c { o\t') . 


— OO 


( 20 ) 
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Due to the Bose symmetry, the two terms in the matrix element (12) give the same 
contribution to the first line of (|20]), and we get 


-?Ht) = 


i d 3 f* K e i ( a+w ')* 
fi 3 (fl + c o' — iA) 


d 3 Kl 


X 


1 

Ft -\- c 


iA £2 + u' 


c f • fci)i0 K ii 2 (fci • ky 

(cuj T oj 1 — iA) 

1 


i d(fj ■ k)*(j)* K ^y g + 


- iA 
iSE, 


){t + 




( 21 ) 


h(Fl + u' — iA) ^2 ' h W ' 0 + u/ —iA' 

where the last line features again a linearly increasing part. This amplitude will be 
used in Sec |3.2| to calculate the average force operator for an atom in constant motion. 
A correction <5c^(t) arising from the acceleration stage is discussed in Sec 
related to the energy stored in the excited state |rf, k). 


4.3 


and 


3. Frictional power and force for constant velocity 

The two-photon power P 2 [Eq.(j9j|] has been introduced and evaluated in detail in 
Ref.@. There, the calculation was performed for a trajectory where the atom is at rest 
for t < 0 and having a constant velocity v for t > 0. We analyze the corresponding 
process in the following Sec|4j and review an alternative calculation reported in 
Ref.® in Sec 5.2 In this section, we evaluate P 2 in the case of constant velocity. 


3.1. One- and two-photon emission 


We find that the only relevant amplitude c 2 2) (t) = ( g , KiK 2 \^(t)) [Eq.(15l] translates 
into the following differential emission rate 

|(c/,KiK 2 |T(t))| 2 d 3 Ki d 3 K 2 


dwo = lim 


£—>•00 


d 4 \ki ■ k 2 \ 2 \f K1 \ 2 \(j) K 


h 4 


2 

1 


7 + {!^ 2 } 


2tt6(uj 1 T co 2 


, d 3 Ki d 3 K 2 


( 22 ) 


n T cuj 

(the factor 1/2 comes again from the Bose symmetry). This yields, to the fourth order 
in d, the power 

d 4 


P> = 


7T 3 h 


I d 3 Ki d\ 2 e ~^ kl+k2)z 


|ki ■ k 2 1 2 
k\k 2 


Q 2 S(u[ + (j 2 ) 

{0 + + w ' 2 ) 2 ' 


x (cai + c o 2 ) Im R{u i) Im R(oj 2 ) 


(23) 
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A simplified evaluation is reviewed in Appendix A.2 leading to a scaling ~ v 4 for a 
small velocity [Eq.(5.4) of Ref.||9|1: 


Po ~ 


512vr 


kv 4 a 2 uj 4 T 2 
8 .,10 


(jJcZ 


(24) 


Such a scaling with velocity was also found within fluctuation electrodynamics BT6l 
H7t f!8l , although the numerical prefactor is different. Eqs.(|23j [24]) coincide exactly 
with one term in Barton's results, called there Pa (see Eq.(5.4) of Ref.f9]|). It is sub¬ 
leading, however, compared to another contribution (called P B ) that scales as Oh 2 ). 
Such a leading velocity dependence was also put forward in Refs.|8j,|34]. We analyze 
the origin of the P B contribution of Ref.0 in Sec j4] where the dependence on the 
atomic trajectory is pointed out. The calculations of Refs.|8i [18] are reviewed in 
Secs |5.2|[53] . respectively. 

With respect to the scaling with the frequency parameters for the material 
response in Eq.(|24]), a similar behaviour has been observed in previous work on the 
metal-vacuum surface where the plasmon resonance is at us = ca p /\[2 


The combination T/uj 2 is then proportional to the specific resistance of the metal. 
Only quasi-DC parameters are relevant for these processes, the spectrum of the 
plasmon pairs being confined to a region of width ~ v/z around zero frequency. 

Since it will play an important role for a generic trajectory [Secj4j, let us also 
discuss here the one-photon power Pi. From Eq.(j8|) we have that to the second order 
in d it is connected with the squared amplitude |cE(f)| 2 / leading to the differential 
excitation rate 


duh = E 


d 2 \ fj • k 
h 2 


-2n5(Fl + u/) d 3 n . 


(25) 


Summing over all final photon states, we recover exactly the excitation rate 


obtained Eqs.(17 181. An exponentially small scaling with velocity still holds 

/Q\ 

when the excitation energy is included in the evaluation of P{ = jdwfhfFt + oj) 
(the superscript indicates again the perturbative order). A consistent perturbative 
comparison with P 2 needs, however, a calculation up to the fourth order in the 
coupling constant. To evaluate the correction P, ( J| to this power in the next order, 
we consider the mixed term 2 Re [cE^CjEf)] and focus on its most divergent part, 
namely the one increasing with t. We find a decrease of the emission rate: 


P, (4) « -7 9 tPi (2) • (26) 

This suggests the resummation P 1 ph P[ 2 ' : ' as expected by the instability of the 
ground state. This shows that, also to the fourth order, the one-photon power is 
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exponentially suppressed leaving, in the case of a constant velocity, P 2 and then a 
force F ~ v 3 as the only relevant contribution to quantum friction. In Sec|I] we 
analyze how these results generalize for a more realistic case where the atom, initially 
at rest, is accelerated to a constant velocity v. 


3.2. Average radiation force 

Before proceeding, it is very instructive to directly evaluate the frictional force given 
in Eq. (j^j) . We consider here the expectation value F(f) = ('E (£) |F(f) |\E(f)) for an atom 
in uniform motion parallel to the surface. We shall use again the expansion of the 
atom+field state |\E(f)) up to the third order of the interaction given in Eq.([To|). The 
nonzero matrix elements of the force operator can be derived from Eqs.( |TT) - [l3| ): one 
just needs to replace the prefactor i in these equations by k or k x 2 . They yield the 
average force in the form 

F(t) = 2Re{ (g, vac|F(£)|? 7 , k)(c* l 1) (£) + c^*(t)^(t) + cf\t)) 


5 ?. 

v 


d 3 nd 3 K 1 d 3 n 2 (ff, K\F(t)\g, (t)^ (t)} , 

where we included products of amplitudes up to order three. 


(27) 


After a straightforward calculation based on Eqs.(14 15 16 21) for the 
amplitudes df' 1 (t) [details in Appendix C |, we find that the average force in the long¬ 


time limit t —>■ oo can be written as F = F (2 ^ + F ( The first term 

TiaQ 


p(2) _ 


7r 


d 3 n k k e~ 2kz Im R{u)6{n + u’) 


(28) 


is second order in the coupling constant and has a simple interpretation: it is the 
recoil due to the emission of a photon. This process is accompanied by the excitation 
of the atom (Cherenkov-Vavilov radiation) and happens at the differential rate du>i 
of Eq.([25]). With every emission, the atom receives a momentum —hk opposite to the 
plasmon momentum. The resulting force acting on the atom is F = — /dwi hk which 


coincides with Eq.(28). 


As explained in Appendix C the fourth-order contribution to the force can be 
presented in the form 

F ( 1J = — -77-P2 — 7 g t F (2) — V v (7 gSE g ) + other exp. small terms (29) 
where 7 g (v) and 8E g {v) are the relaxation rate and Lamb shift of the ground state. 


respectively, Eq.(17). Recall that "f g (y) arises from quantum Cherenkov-Vavilov 
radiation and is exponentially small. This is also true for the second-order force F (2) 
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[Eq.(|28]>] because the resonance condition 9 + oj' = 0 involves the same threshold as 
Eq.([T9]) for The 'other exp. small' terms not written explicitly in Eq.(|29]) have a 
similar origin. 

The important result is that the force (29) gives the leading order for small 
velocities and involves the two-photon power P 2 obtained in Eq.(|23]). This means 
that for a uniformly moving atom, the average radiation force starts like 0(v 3 ), 
which is coherent with the radiated power obtained in the previous section. The 
force calculation thus provides an independent confirmation that two-photon rather 
than one-photon emission (plus atomic excitation) is the dominant loss process. A 
comparison with the results of Dedkov and Kyasov 113911 and Intravaia et al. Ifl8ll is 


made in Appendix A.2 agreement up to a numerical factor is found when the atomic 
transition is off-detuned with respect to the surface plasmon resonance, O C w s . The 
dependence on distance involves the steep power law F ~ 1 /z 10 . 


4. Accelerating the atom and subsequent radiation 

In this section, we consider atomic trajectories that are accelerated over a finite 
duration before reaching their final velocity v. This material generalizes the 
calculation of Ref.0 of the two-photon process where a second term (called P B ) in 
the two-photon emission was found that scales with 0(v 2 ) in velocity. The main 
result is that the term P B depends sensitively on the way the atom is accelerated. 
Ref.|9j is only recovered for a 'sudden boost' (infinitely short duration), while in 
the opposite or 'adiabatic' limit, P B becomes strongly suppressed. The scaling with 
velocity 0(v 2 ) is maintained, though. 

To interpret this behaviour, we have also evaluated the one-photon power P\ 
and found, quite surprisingly, that for accelerated trajectories it is not exponentially 
small (as in the previous section), it is negative, and exactly cancels the two-photon 
emission P B . This suggests the following picture: the acceleration stage creates a 
finite occupation p e ~ v 2 of the excited state (including one photon). The excitation 
process is qualitatively similar to the 'acceleration-induced radiative excitation of 
ground-state atoms' analyzed by Barton and Calogeracos |38|. Subsequently, this 
'real' rather than 'virtual' excitation decays resonantly into another photon. The 
resonance condition fixes the energy of the second photon, so that the radiative 
power captured by the term P B ~ p e ^ e hFl is balanced by a decaying excitation 
probability (negative Pi). 

The calculation proceeds by working out the probability amplitudes, starting 
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(t)= ~- 


d(v ■ k*)(j>* 


dti e 


i(0+aj)4i c -ik r(ti) 
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(30) 


where the matrix element ( [TT] > of the atom-field coupling was used. Recall the 
compact notation n = kca for the photonic modes and note that we have kept a 
generic atomic path r(/i) under the integral. The / |-integral appearing here will be 


denoted A(e, k: t ) and discussed in detail in Sec, 4.1 We prove there that at large times 
(once the launch is completed), the amplitude takes the form 


t » r : Ci^(f) « - 


d(ff ■ k*)f* K f e 1 ( n+a, 0* 


+ B, 


e,K 


(31) 


h L i(0 -h co') 

where the first term is the same as for a constant-velocity path [Eq.(|l4))]. We 
interpret the second term as a non-adiabatic excitation process whose amplitude is 
approximately (for small velocity) given by 


i fk • v) e 'k-r(o) 

B e ,. ~ ^ ^ S((0 + 0J)T) . 


( n + u )» - (32) 

The dimensionless factor E((12 + u)r) depends on the specific shape of the path. It is 
proportional to the Fourier transform of the acceleration [Eq.(|50])] and decays to zero 
when the product of duration r and frequencies 12 + u> is much larger than unity. In 
the opposite limit ('sudden acceleration'), E((12 + u)r) —>■ 1. 

In the next order of perturbation theory, we deal with the two-photon amplitude 


L 

4 2) (t) = j dt 2Ci ) {t2){g-,K 1 ,K2\V(t 2 )\f};K) 


(33) 


= ^(k, ■ £,)*«, / di 2 A(e, Kl ; t 2 ) + {1 « 2} . 

—oo 

The additional term denoted {1 <(->■ 2} makes this expression symmetric under 
plasmon exchange. The t 2 -integral written here will be called Ai\ 2 (f) similar to 
Ref.|9l. We find (Sec|4.2[) for this two-photon amplitude the asymptotic form 


t » r : .AT 12 (f) 


3 i(a)^+0)2 )t 


+ V 


B, 


e,K 1 


+ cst _ ; (34) 


(11 T ca()(ca( T ca 2 ) i(ca 2 — 12) 

where the cst. denotes /-independent terms. The first term again recovers the 
previous constant-velocity result from Eq.([l5]>. The second term is proportional to 


the non-adiabatic excitation amplitude £> ejK [Eqs.(31 32)] that appeared in the first 
order. 
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The power of two-photon emission from Eq.Q is proportional to \M\ 2 it) + 
.M21(f) 1 2 /t in the large-t limit. This is calculated in Eq.([45| below. In Sec j4.3[ we 
finally discuss the rate of change of the energy stored in the excited state \ff, k) (the 
one-photon power Pi defined in Eq.Q) and show that it balances exactly the 0(v 2 ) 
contribution to the two-photon power. 


4.1. Exciting the atom: the one-photon amplitude 


Let us consider the first step of the physical process described above. The one-photon 
transition amplitude is proportional to 

t 

A(e, k; t) = fdtie^e-^. (35) 

— OO 


A few general properties of the first-order amplitude can be secured without 
specifying a particular path. We only require that r(t) changes its velocity around 
t — 0 with a typical duration r. We also assume that the origin of the coordinate 
system is chosen such that for i>r ; we have r(£) ~ vf (see Fig|3]for a sketch). Let us 
focus first on t > 0. We split the integral into — 00 < t\ < 0 and 0 < A < t r leading to 
a natural decomposition A(t) = A-+A+(t). For the term A+(t), we perform a partial 
integration after subtracting and adding vfi in the exponent [see Eq.(|30j)]. This leads 
to 


A + (t) = — 


-,i(Q+w)tl ik-r(ti) 


i(0 + to — k • v) 


+ / dti 


k • (r — v) 

+ uj — k • v 


_,i(f2+oj)ti —ikr(ti) 


(36) 


The advantage of this representation is that the first term yields what we call the 
'adiabatic limitj|] 

gi (Si+u)t g-ik-r(t) 


k) t) — - 


(37) 


i(0 + ui — k • v) 

which is nothing but the term we found previously for a trajectory with constant 
velocity [see Eq.(|l4])]. In the integrand in Eq.(|36|), the difference r — v vanishes as 
soon as the atom has reached a constant velocity. Hence the integral approaches a 
constant for t r. In this limit, we therefore proved that the excitation amplitude 
takes the form 


t»r: A(e, k; t ) = A+ d) (e, k; t ) + B C)K + o(l) , 


(38) 


I This is not a fully adiabatic expression since the denominator contains v instead of the instantaneous 
atomic velocity r (f). 
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Figure 3. (Color online) Three possible atomic trajectories: smooth continuous boost 
(thick black line), linear acceleration ramp (thin blue line), and instantaneous boost 
(dashed gray line). See also the table below for more detail. 
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x{t) 
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0 for t < — r 

(t +t)w/(2t) for — r < t < r 
v for t > t 

jo for t < 0 
] v for t > 0 


0 for t < —r 

v/{2 t) for — r < t < r 
0 for t > t 

vS(t) 


where B CyK can be read off by adding A_ to the remaining terms in Eq.([36]). The error 
o(l) in Eq.([38| is made of terms that vanish in the limit i>r. A similar manipulation 
can be applied when t < 0. No splitting and subtraction are needed, and we get 

i(Q+uj)t -ik-r(i) t k . f 

t< 0: A(e,K]t) = 


i(tl T ta) 


dfi^^e l(J]+u)il e" ik ' r(tl) .(39) 

+ U) 


The first term may again be called 'adiabatic' and does not contain any contribution 
from the lower limit t — —oo because we assume that the atom-field coupling 
is switched off there. Note also that there is no Doppler shift in the frequency 
denominator. The second term vanishes for t <C — r when the atom is still at rest. 


4.2. Emitting the second photon 

Here, we focus on the two-photon amplitude in the long-time limit t —f oo. Using 
the time scale r for the 'acceleration stage' of the atomic trajectory, we assume more 
specifically t>r. Let us introduce a time t a with the property r « « i such that 

the atomic velocity is r(£ a ) = v. We split the integration range of the C-integral in 
Eq.([33|) into —oo < t 2 < t a and into t a <t 2 <t obtaining Mi 2 (t) = M- + M + (t) (and 
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similarly for M. 2 i(t)). We shall see that only the probability |Al + (t)| 2 contains terms 
growing linearly with t and then contributing to the radiated power, while the rest 
tends towards a constant. The contribution M + {t) can be evaluated by inserting the 
asymptotic form ( |38| > for the amplitude A(e, K \; f) into the integral: 

t 

M+(t) = I d t 2 k i5 t 2 ) + B eiK1 ) 

t a 

= M { +\t) + M { +\t). (40) 


Because, for t > t a , the atomic velocity is constant and the first term evaluates to 

pi(‘d+‘*4)i2 t 


M?\t) = 


(w' = uJi - ki • v) 


(41) 


(12 + Sj'fjibj'i + Oj'f) 

we recover the result for a constant velocity [cf. Eq.(|l5])]. We symmetrize under 
photon exchange, square and get in the limit t t a : 

(211 + + uj' 2) 2 


M^it) + {1 -H- 2} = 2n(t - t a )+ u'f)- 


(42) 


(12 + o >() 2 (12 + u 2 ) 2 ’ 

where the sum + co 2 in the numerator can of course be omitted. The second term 
in Eq.(40l is an elementary integral as well: 




(t) = 


B, 


' e .«l i(w'-f i)t 2 


\{u' 2 - 12) 

We symmetrize again and identify the individual squares as the leading terms: 

\t) + {1^2} = 2n(t — t a )5(uj 2 — 12)|23 eiKl | 2 + {1 -f-)- 2} . 


(43) 


(44) 


It can be checked that the 'mixed terms' in the squared amplitude lead to 
contributions that either oscillate or tend to constants as t —> 00 (no inverse squares 
like 1 / {'jj\ + uj' 2 ) 2 appear). This also holds true for the mixed terms involving one 
factor M.- (see Appendix B for more details). Finally we get 


lim 

t— i>00 


\-M-12 (t) + M 2l {t)Y 


= 2'k8(u' 1 + u 2 ) 


412 2 


(12 + ca() 2 (12 + ui'f) 2 


+ 2nd(uj 2 — 12) \B e ^ Kl f + {1 -<->■ 2} 


(45) 


From the above calculation, one can see that the first line of Eq. ([45]) involves only 
the constant-velocity part of the atomic path and hence does not depend on the way 
the atom is put into motion. As was reviewed in Secj3]and Appendix A.2 this term 
contributes the amount Pa ~ v 4 to the two-photon emission. The second line of 
Eq.([45|) leads to 

■ 2 12 2 '■ 


Pb = 


a 


Ah 


d d K 2 (wi + cu 2 )|/ci ■ k 2 \ 2 \f 1 (j) 2 \-2n8(uj 2 - 12)|H e> « 


(46) 
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and scales with v 2 for small velocity [Appendix A.l[. This arises because the 
amplitude B e , K ~ v, as we now show. 


We start with the general expression collected from the terms in Eqs.(36 39) that 
become constants for large t > r: 


B = 


e -ikr(0) 

i(14 T ca) 


o 


df | ^ ' ** c i(n+o;)U g —ik r(ti) 


Ft + UJ 


a - ik-r(O) 


i(14 + u — k • v) 

The first terms in the two lines sum to 

e -ik-r(0) g-ikr(O) 


f dti J k '^ ^ e m+w)tl e~ ik ' r(tl) . (47) 

14 + uj — k • v 


-k • v) e lk r (°) 


i(14 + o;) i(14 + uj — k • v) i(14 + ca)(14 + uj — k • v) 

This is also what is found for Ref.0's 'kink trajectory' where r = 0 for t < 0 and 
r = v for t > 0. The remaining integrals are bounded by 0(kvr/(Fl + u)) so that 
their contribution to the amplitude £> e , K is also linear in the velocity v. To proceed, we 


apply another partial integration to the two integrals in Eq.(47). Summing the results 
and expanding for small v, one gets 


k • r 


k ■ (r — v) 


" ~ c i(n+ 0 ))ti g —ik-r(ti) _|_ J ^ 11 V > gi(0+^)ti g-ik-r(ti) 


Ft A oj 


k • v e lk r (°) 


i(fi 


UJ 


g-ikr(O) 

i(fl + co) 2 


'FI + uj — k • v 
u 

df 1 k-fe i(f2+aj) * 1 


(49) 


where the first term cancels with ( [48] ) to leading order in v. We are thus left with the 
Fourier integral of the atomic acceleration 


B f , 


ikr(O) 


i(fl + uj) 2 


CXJ 

f dfik-r(f 1 ) e i(f2+aj)il = 


i (k ■ v) e 


—ik-r(O) 


(14 + uj) 2 


-S((14+cu)r) .(50) 


which is the result announced in Eq.(32). This already permits us to draw a 
conclusion for a generic trajectory whose velocity is monotonously raised to its final 
value v. The Fourier transform of the acceleration then exists and is maximal for 
14+ca = 0. If the acceleration occurs over a finite duration 0(t), the Fourier transform 
drops to zero when the frequencies satisfy (14 + uj)t 1. On physical grounds, it 
seems quite plausible that the acceleration stage takes more than a few femtoseconds 
while 14 is typically in the visible range (and uj > 0). This inequality is therefore 
amply satisfied, and the corresponding two-photon emission is strongly suppressed. 
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It is interesting to examine different atomic trajectories r (t) and provide a 
quantitative analysis how the 'acceleration stage' affects the final result for the two- 
photon emission P B ■ A trivial example to begin with is an inertial path with constant 
velocity. The acceleration is zero at all times, and Eq.([50)) gives B,,, K = 0. Hence the 
two-photon power P B = 0, and only the 0(v 4 ) contribution called Pa from Secj3] 
remains. 

Our second example is Refill's instantaneous boost, r(f) = v5(t). The trajectory 
is plotted as a dashed line in Fig[3j The Fourier integral gives 

I B I 2 « ~ v ) (51) 

\ u e,K\ /o i 

(iz + cu ) 4 

This result can also been inferred from Eq.(4.10) of Ref.0 by writing it in the form 
given in Eq. ([45]) (and taking into account the erratum). 

As a third example, consider a linear velocity ramp, as illustrated by the middle 
path in Figj3} Velocity and acceleration are given in the caption: specifically the 
acceleration is constant during an interval of length 2 r (details of the full calculation 
for this path can be found in Appendix Bj ). The sinc-function resulting from the 
Fourier integral (|50|) then gives 


I B e . 


(k-v) 2 sin 2 [(fl + ui)t\ 


(52) 


(fl + cu ) 4 (fl + cu) 2 r 2 

The second fraction reproduces Ref.0's path in the limit r —y 0, but gives a strong 
reduction in the opposite case. Since the typical frequencies contributing to the 
integral (46) are u)\ ~ lus due to the plasmon pole, \fi\ 2 ~ Imli(wi), the power P B 


gets reduced by a factor l/(r(fl + co s )) 2 -C 1. This result for the linear velocity ramp 
has been reproduced through a differently routed calculation by G. Barton (private 
communication). 

Finally, let us consider the 'smooth boost' plotted as a thick black line in Fig[3] 
The acceleration is given by an infinitely differentiable function (see figure caption) 
whose maximum value is v/(4r) and whose width is O(r). Evaluating its Fourier 
transform, we get the probability 

-i 2 


| b I 2 = 


(k-v) 


7r(H + lu)t 


(53) 


(O + ta) 4 |_sinh[7r(n + u)r} 

The second fraction in this expression shows that compared to the 'kink path', the 
power P B becomes exponentially small in the limit of an adiabatic boost rQ 1. 

Fet us attach a physical meaning to the quantities calculated in the previous 
sections. Going back to the Schrodinger picture, the first-order amplitude c^\t) 
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(54) 


The first term, independent of t, can be understood as being part of the ('dressed') 
ground state \q, vac) (still at zero energy), where the atom is surrounded by a 
('virtual') cloud of photons (plasmons). The second term oscillates at the (bare) 
energy of the excited state \ff, k), but including the Doppler shift (a/ instead of uj). 
It can be shown that the Hamiltonian of our atom+field theory can be transformed 
to a time-independent form by going into a frame moving with the atom. (Details 
are postponed to another paper.) In this picture, the state | ff, k) evolves freely 
at the frequency Fl + to'. We therefore conjecture that along a path with a time- 
dependent acceleration, the amplitude describes the 'real' excitation of the 
atom+field system | |38| . The required energy transfer is in heuristic agreement with 
the frequency uncertainty arising from the finite duration of the acceleration, as 
expressed in the Fourier integral ( [50] ). 

We can also define an excitation probability (not a rate), summing over the 
plasmon states and the three sublevels | if) 



v 



(55) 


where the ^-integral can be performed, yielding 3nv 2 / (4,~ 5 ). 

4.3. Excitation power: a subtle cancellation to the fourth order 

The previous analysis examined in detail all the components of the physical process 
describing the acceleration, the excitation and the subsequent radiation of an atom in 
motion near a surface. This was necessary in order to discern and generalize, to the 
case of a generic trajectory, each single contribution to quantum friction. In Ref.|9j the 
friction force F is calculated based on the identification with the two-photon power 
loss, P 2 = —v • F. For a generic trajectory this calculation led to P 2 = Pa + Pb 
where Pb ~ 0(v 2 ) gives the leading order. Barton thus concludes that F ~ v for 
small velocity 0, at least for the trajectory called instantaneous boost above. In 
the last section we showed, however, that the value of Pb depends on the detail of 
the trajectory: a smooth boost gives a significant reduction, and a constant velocity 
simply leads to Pb = 0 (see also Secj3]). 
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Although physically sounding, the calculation based on P 2 is incomplete since it 
does not take into account the power needed to create the excited state \ff, k), called 
the one-photon power Pi in Eq.([8]). A similar omission in earlier work was criticized 
by Volokitin and Persson, see Ref.|[T3f. (For an analysis of the 'internal energy' 
appearing in Pi, see also Ref. lf39i for example.) The sum Pi + P 2 translates the change 
in the total energy of the evolving state |T(£)). This energy is not conserved, since 
the interaction is time-dependent. For the same reason, only the free Hamiltonian 
H 0 = Ha + H F (atom and field) is used to define the energy of the state |T(£)): 

(*(t)\H 0 \*{t))Kt(P 1 + P 2 ) + ... (56) 


The power P\ is calculated again by pushing to the third order the one-photon 
amplitude (if, ac|\ 1/(£)), called c^(f) in Eq.( 10). This extension is necessary because at 


first order, even with an acceleration phase, the excitation rate is exponentially small 
[see Eq.(|l9|)]. We focus again on the state sequence |p, vac) —> \if,n) —> \g, KiK 2 ) —> 
| ff, k) passing via the two-photon state. (The sequence via the ground state \g, vac) 
gives again exponentially small contributions.) 

The calculation proceeds along lines similar to Secj3j Perturbation theory yields 
the integral (|20|> where we insert now the expression for (g, /ci« 2 |\l/(f)) generalized to 


the case of a generic trajectory (proportional to the amplitude M(t) from Sec 4.2). In 
addition to the constant velocity result we obtain a correction (t ) to the amplitude 
coming from the second term in Eq.(j34) and its symmetrized partner. We consider 
that interaction times in the interval r < f 3 < t give the main contribution and 
approximate 


d 3 

h 3 


tic ( i\t) PS ^3 ld 3 Ki(ki ■ ky\f K1 \' 2 (f)* K (ff■ ki) 






+ 


B P 


(57) 


i(o/ — 0 — iA) \{u'i — Ft — iA) 

The second term in the curly brackets gives rise to a linear increase in time. The 
fourth-order approximation to the excited-state probability | (if, k\ T(£)) | 2 ps |c^ (t) \ 2 + 
2 Re[ci 1) *(£)5ci 3) (£)] + ... thus provides us with an excitation rate 

i® K |<i<(«))i 2 -i4 1> mi 2 


lim 

t—> OO 


2<* 4 „ 
“^ Re 


(if - k)\(j) K \ 2 / d 3 ni(ki ■ k)*\(J) Kl \ 2 ( r if ■ kf- 


I B I 2 
I | 


(58) 


i(ca / 1 — Q — iA). 

Note that for this rate of change, we only need the non-adiabatic amplitude B e>K in 
the first-order expression c 1 , 1 j *(£) [Eq.(31|], since the other combinations give rise to 
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oscillating contributions. The one-photon power (|8j) becomes 


Pi = 


a 


Ah 


d 3 Kd 3 Ki(Fl + u)\(f) K \ 2 \(j) Kl | J |• k\ 2 \B e , K \ 2 2'r[5[uj' l — Ft) . (59) 


To the lowest order in velocity (recall that B r , K is proportional to v), we may drop the 
Doppler shift in 5(ca( — Q). Now, it is possible to check that the one-photon power (591 
exactly balances the contribution P B to the two-photon power [Eq.(46)]. It leaves Pa 
as the only relevant contribution for the total dissipated power, even for a particle 
path including an acceleration phase. 

The work developed in this section is the central result of our paper. Firstly, it 
shows that the perturbative approach described in Ref.U9J strongly depends on the 
acceleration phase that brings the atom to a constant velocity v. Secondly, it proves 
that the description given in Ref.O of the quantum friction is incomplete and that, 
when corrected, it is in agreement with a drag force at zero temperature proportional 
to v 3 . 


5. Results from fluctuation electrodynamics 

In the previous sections we provided a complete description of quantum friction 
within the framework of perturbation theory. This approach has the merit of 
relying on well-established techniques, even if the mathematical machinery is 
somewhat cumbersome. Quantum friction, however, has been examined within 
other frameworks, approaching the problem from other perspectives. For the sake 
of completeness, we review in this section some of the results from fluctuation 
electrodynamics, which is one of the most used approaches to describe the quantum 
mechanical interaction of two neutral objects. 

5.1. Spectral densities 

Correlation functions of the atom and field variables are a convenient way to 
characterize the atom-field interaction in terms of 'resonant' and 'non-resonant' 
processes. We start by collecting a few formulas for the free observables and evaluate 
their correlations in the 'bare' ground state denoted as \g, vac). 

Field correlations. We use 'time-ordered' correlations as they often appear in time- 
dependent perturbation theory. For the free scalar potential and t > t' 

C$(r, rt — t') = (vac|l>(f, f) < l(r / , t')|vac) 
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pj rJc r 

= w I IT J iul ImBMe 

0 


ik(r-r')-k(z+z') e ~ioj(t-t') 


(60) 


Evaluating this for an atom with constant velocity v, we get r(t) — r{t') — (t — t')v 
and observe that the vacuum spectrum extends into negative frequencies, of the 
order 0 (—k ■ v) = 0[—v/z). This estimate is based on the natural momentum 
cutoff provided by the exponential e~ k ^ z+z '\ The rest of the frequency dependence is 
governed by the reflection amplitude R(uj ): a peak at the surface plasmon resonance 
so = cos with width T and an algebraic decay ~ 1/ta 3 in the UV. In the time domain, 
these features translate into a correlation that oscillates at ujs with an exponential 
envelope of width 1/T, plus an algebraic long-term tail ~ l/(t — t') 2 that arises from 
the 'Ohmic' behavior Im R(oj) ~ u for cu —> 0. For the electric field, evaluated along 
an atomic path parallel to the surface, we get similarly (u' = u — k • v, frequency in 
the co-moving frame) 

(vac|£'t(r(f), t)Ej(r(t'), f 7 )|vac) 

h 

~ 2^ 


U.U* P ~ 2kz °? 

d 2 k 1 j - - [ du/ Im R[u' 


k 


k • v)e 


—i 


(61) 


-k v 


Note that for a more general trajectory, the correlations are not stationary, and 
more involved spectral representations like Wigner or wavelet transforms would 
be needed. The response function of the free field is known as the Green function 
(tensor). Standard linear response theory gives 

Gij(r,r',t-t') = ^<d(t — t 1 ) (vac | [Ei (r, t), Ej (f\ t ; )] | vac) 




-k(z+z') 

®_pik(r-r') 

k 


x Jdu Im R(u) [kik* - kjk* 


(62) 


with an obvious evaluation along the path of a moving atom. (It can be checked that 
the last line of Eq.([62)> agrees with the solution of the Maxwell equations for a point 
dipole in the non-retarded approximation.) 


Dipole correlations. The free dipole operator shows in the theory of Ref.[[9] a sharp 
line. It is actually the specific challenge of this model that the line broadening appears 
self-consistently at higher orders in perturbation theory. In the atomic ground state 


<g\D,(t)bM)\g)W = 6, 


haD 


u 


s —i 


2 


(63) 
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after summing over the degenerate excited states \ff). In a simple scheme where the 
states \fj) have decay rates 7 fj, this could be generalized to 

haFl 


{g\D t (t)D s (f)\ 9 ) ■■ 

giving a Lorentzian spectrum: 


E ViVj 




EH = / dre lu)T (g\Di(t + r)Dj{t)\g) = —— E 


IvViVj 


(u - FI) 2 + 7 S /4 ' 


(64) 


(65) 


This is also the result of master equation techniques in combination with the 
regression formula Il22l . Fermi's Golden Rule yields for the decay rates in front 
of a smooth metallic surface [Eq.(2.10, 2.11) in Ref.lf9lll: 7 7 = rgq^rjj 7 , 7 = 
(af2)/(4z 3 ) Imf?(fl), where q,j is a dimensionless diagonal tensor with elements 
q xx = q yy = 1/2 and q zz = 1 . We recognize again that Imt?(o;) gives the spectral 
density of the plasmon field. 

The atomic response is given by the polarizability tensor whose lowest 
approximation in the spectral domain is 


«i 0) H = i /d re iulT {g\[D i (t + T),D j (t)]\g) 


h 

0 

_ c 

5ij n 


+ 


SijaFl 2 


f2 2 — (u + iA) 


2 > 


( 66 ) 


^Fl — oj — iA Fl T ui T iA, 

where A is a positive infinitesimal that shifts the frequency into the upper half¬ 
plane. This has the same structure as for an oscillator. A simple finite-damping 
generalization would replace A by the actual line widths: 

ViVjOiFl 2 


qgM ~ E 


Fl 2 


Ua 


+ i7 ffY 


(67) 


(For a discussion of the imaginary part near the anti-resonant peak fl + u « 0, 
see Refs. li23l l24il d We note that for infinitely narrow lines, the dipole correlation 
functions in Eqs. (63 66 ) satisfy the zero-temperature fluctuation-dissipation (FD) 
relation 


Sij(co ) = 2hQ(u) Im Oiij(u) , (68) 

where S l j{uj) is the dipole correlation spectrum in Eq.([65]). The FD relation is not 
satisfied, however, by the line-broadened expressions presented in Eqs. ( 64|67 l: the 
dipole spectrum does not vanish like Im (ca) ~ cu near zero frequency, and extends 
also into the negative frequency band. In general, however, these expressions are the 
result of approximations. We recall that the FD relation is valid under relatively mild 
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equilibrium requirements, in particular it also holds when the dynamics of the dipole 
operator is non-linear H26ll27ll . For driven systems like in our case, generalizations of 
the FD relation in Eq.([68]) |28f involve additional 'source' terms |1T8| or correlations of 
observables that are conjugate with respect to entropy (production) rather than the 
Flamiltonian 11291 . 

5.2. Macroscopic QED with Markov approximation 

Scheel and Buhmann derived the quantum friction force on an atom of arbitrary 
internal state from the average Lorentz force 0, which, in the non-retarded limit, is 
given in Eq. (|3j) at the beginning of this paper. The time evolution of the electric field 
is obtained by formally integrating the equation of motion for the bosonic operators 
a K , of with the result that E(r, t) = E iree (r, t ) + E lS) (r, t) with the free field operator 

OO 

J free (f, t) = fduj £(r, u) e~ iujt + h.c. . (69) 

o 

When the atom is not externally driven, then we find in normal ordering and an 
initial vacuum state for the field, (• • -S(r,u)) = (<D (r, ui) ■ • •) = 0. When evaluating 
the radiative force F(t) in normal order, it thus turns out that it is entirely due to 
radiation reaction, i.e., the relevant electric field is the source field emitted by the 
atom at previous times. This can be written with the field's Green function [Eq.([62]>]: 

E\ S) (x, t ) = j df' Gij(x, r{t '), t - t')bj(t') , (70) 

where D{t') is the dipole operator. The velocity-dependent force is due to the delay in 
the radiation reaction field: the atom acts as a source for the electric field at a previous 
point on its trajectory r(t'); the generated field then causes a force at a later position 
r(t) where the atom has moved to in the meantime. For an atom moving normal to 
the surface, Doppler shifts of atomic transition frequencies and line widths give rise 
to additional velocity-dependent effects. At retarded distances, the Rontgen coupling 
of the moving atom to the electromagnetic field needs to be taken into account 0. 

For our problem with short (non-retarded) distances to the surface, the Green 
function in Eq.([62j) yields a natural split of the source field into positive and negative 
frequency components M s \x, t) = £ iS) (x, t) + h.c., where 

-* h °° 

£ (s) (x,t) = ^ fd 2 kye- k{x3+z{t '» [dw Im R(u) 

2tt 1 J k J 

o 

t 

J dt'e ik ' (x - r(t,)) - iaj(i - f,) (F ■ D(t')) . 
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—OO 


(71) 
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In the following, we evaluate this at the position x = f(t ) of the atom and assume that 
the latter is moving at constant velocity v parallel to the surface. Although this is not 
the most general trajectory, we will later argue that within the Markov approximation 
used in this subsection the precise history of how the particle achieves its terminal 
velocity does not matter. The average Lorentz force thus becomes 

F(t) = (A(i)V£f , (> s '(«),«)+h.c.) 

1 h °° 

= - 2?/ d V ; d “ Imfl( " ) 

0 

t 

x J +c.c. . (72) 

— OO 

In Ref.f8], this expression was expanded for small v; for the ease of comparison with 
other approaches, we defer this to a later stage [Eq.([77j) below]. 

For weak atom-field coupling, we may evaluate the dipole-dipole correlation 
function (Di(t)Dj(t')) using the Markov approximation. This entails converting the 
equations of motion for the atomic operators into an integral equation and taking 
the slowly-varying operators out of the integral. The result is an effective solution 
to the equations of motion involving only operators at a single time t. Hence, all 
memories of previous quantum states have been lost. As shown in Ref. |40| , the 
Markov approximation may become invalid, e.g., if an excited atom near-resonantly 
interacts with a narrow resonance of the medium-assisted field (which is not the case 
here). For our case of an atom initially prepared in its ground state, the upshot of this 
analysis is the following representation of the dipole correlation function in terms of 
lowering operators A gr j between the atomic levels: 

(g\Di(t + T)bj{t)\g) = d 2 J2ViVj(9\Agfj(t + T)Ab(t)\g) (73) 

v 

-^{g\A gi f(t + r)A ] gi .(t)\g) = ( - i ft- ^)(c/|i 9 , 7 (t + T)A ] gfj {t)\g) , (74) 

where the second line contains the atomic frequency Ft and the line width of the 
| g) if) transition. This also yields the correlation function given in Eq.([64|>, using 

in the initial condition the closure relation ^^(glAg^fAbft^g) = 1 for a ground- 
state atom. Coming back to the radiation reaction force [Eq. ([72])], at large times the 
integral over t! evaluates to 
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J dt' e <p IA (0 A (*') I <?> =<?£• 
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(75) 
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F = 


kafl 

(2t r) 2 


r d 2 k^e~ 2kz 


dca Im -R(ca) 


\v-k\ 2 ln 


(Q + a/) 2 + 7^/4 ’ 


(76) 


0 V ' ' '»?/ 

where a/ = ca — k • v is again the Doppler-shifted frequency. If this were evaluated for 
an infinitesimal (and isotropic) linewidth, we would recover the first-order force from 
Eq.([28]), exponentially suppressed for small v. Following Ref.El, we keep a finite 
linewidth, observe that the lateral force vanishes for an atom at rest, and expand for 
small v : 


P ~ _ 


hafl 
~ 2 n* 


d 2 k 


k(k 


—2kz 


k 


Y \n- k\ 2 lfj /dca 


Im f?(ca) 
(fl + ca) : 


(77) 


v o 

We recall that this result holds for an atom moving parallel to a surface at nonretarded 
distances. Note that due to the Markov approximation made, no memory of previous 
times is retained in the evolution equation for the atomic variables. In particular, 
this means that this result does not depend how the atom is accelerated to its final 
velocity v. The line widths ~<fj for the smooth metal surface of the present model 
have been given in the previous section where their anisotropy was also discussed. 
Note that this and the friction force were incorrectly given in Ref. 0 due to an 
error in the averaging over excited states. The corrected calculation can be found 
in Ref. ||4Tfl. The sum over the excited states, weighted with their line widths, leads 
to Yjfj\v ' kflrj = lyfc 2 , where 7 = is the line width parameter for a perpendicular 
dipole. The frequency integral in Eq.([77[) can be performed with a Wick rotation to 
the imaginary axis 


dca 


Im R(lo) 
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7 Q 2 _ 2 
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(78) 


0 v .ca) 3 J "(E ! 2 + a 3 4ca 5 (T! + u; 5 ) 3 

where the last expression was obtained for a narrow surface plasmon resonance 
(T os, see also Eq.( |A.2[ )). Performing the /. -integral, we finally get for the lateral 
force 

3 UaFl 


F = — 


Upl 


I67r cus(f2 + cus) 3 
9 ha 2 Fl 3 


Jd 2 kkk(k ■ v)e 


—2 kz 


= — V 


S 4r 


512^ 8 us{n + ca s ) 3 (^ 2 - ca|) 2 ’ (79) 

which gives a frictional power — F • v which agrees with the value for P B given 
in Eq.(|46]) and first derived in Ref.[[9l. The significance of this agreement remains 
unclear at the moment due to the very different underlying assumptions. The 
calculation reviewed in this subsection depends only on the final atomic velocity. 
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the details of its launching procedure being lost in the memory-less Markovian 
behaviour due to the finite correlation time resulting from atomic dissipation 
(spontaneous decay). On the other hand, the time-dependent perturbation theory 
of Ref. K9J is valid for small times and so implicitly assumes an infinite correlation 
time. It hence depends on the atomic acceleration trajectory, where the agreement 
with the above result is found only for a very specific out of many possible choices: 
sudden acceleration. For a more meaningful comparison, a calculation along the lines 
of Ref. 0 could be performed for a dissipative system with a finite correlation time, 
where at sufficiently large times one would expect the result to also be independent 
of the acceleration stage. 

5.3. Non-equilibrium dipole correlations 

The approach followed by Intravaia, Behunin and Dalvit Ifl8l combines techniques of 
fluctuation and macroscopic electrodynamics. While the expression for the radiation 
force has the same structure as Eq.([72]) above, the calculation of the dipole correlation 
function is performed differently. In the limit t —>• oo, the system becomes stationary 
and the correlation function 

Cij(t, t — t) = (Di(t)Dj(t - t)) -> Cij(r] v) (80) 

depends only on the time difference t — t — t' and the final velocity v. (Corrections 
due to the acceleration stage drop out at this point.) In the previous expression the 
dipole operator D(t) contains the exact dynamics of the moving atomic dipole (all 
orders in perturbation theory), i.e. including the backaction from the field/matter. 
The mean value has to be evaluated with respect to p NE ss = lini/^oc p(f) which defines 
the (in general unknown) density matrix describing the non-equilibrium steady state 
(NESS). The latter obviously depends on the atom's velocity v; which is why we 
added the second argument v to the correlation function. 

For a dipole operator with the structure D(t) = d J2ffV{A g f}(t)+A' gf .(t)), the matrix 
Cij (r; v) is symmetric, and since stationarity implies Cij (r; v) = Cf (—r; v), the power 
spectrum 

/ OO 

dr e ioJT C , u (r: v) (81) 

-OO 

is symmetric and real. The frictional force can then be written as 

1 \c 00 

F = — —- J d 2 fc — e~ 2kz j dca Im R(u)kik*Sij (k • v — ce; v) . (82) 
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In order to evaluate the previous expression one needs to know S t j(uj\ v), which is in 
general available only within a perturbative approach. (An exception is an isotropic 
oscillator atom for which the dipole power spectrum can be found exactly |iT8|.l 
For the model atom of Figjljright), there is a Pauli algebra for each excited state 
| if), spanned by the operator = A g g + A^ gfj = \rf)(g\ + \g)(rf\, together with 
o 2 = i(\g)(ff\ — \ff)(g\) and <r 3 = \ff)(fj\ — \g)(g\. With the atom+field coupling 
Vit) = *-Di(t)Ei(r(t),t), we have the following nonlinear equation of motion in the 
Heisenberg picture 


+ f2 2 <7i(t) = 


2dn 


a 3 (t) ff- E(r(t),t). 


(83) 


We focus our attention on the computation of the two-time correlation tensor 
Cij(t,t';v ) = d 2 v ) with C(t, t'\ v) = (<7i(f)<Ti (£')). To the lowest order 

in d, it can be evaluated from the free evolution of the dipole operator, resulting in 
C-j\t , £'; v) = d 2 5ij e -' n (i-A which is nothing but Eq.(63>. This results, however, in a 
frictional force that is exponentially suppressed in 1/v (see also the discussion after 
Eq.(|76|)). To get a force scaling as a power law in v, one needs to include second-order 
radiative corrections in To this end we first insert in Eq.([83]> the formal 

solution for the dynamics of 


2 dr- ~ 

d 3 (t) = <7 3 (—oo) + — J dtid l (t 1 )ff-E(f(t 1 ),ti), 


(84) 


and then replace the exact field E(r, t ) by its free evolution E free (r. t ), given in Eq.( 691. 
This leads to an equation of motion correct to the second order in atom-field coupling: 

*2d 2 r * r A , 

(*) + -rr / dt lVi r,j{Et™(r(t),t), Ef ee (r(ti), 
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(85) 


Multiplying this equation from the right by aft 1 ) and taking the expectation value 
over the initial state |g,vac) (we recall that the bare initial state can be used here 
because corrections to the NESS are captured in perturbation theory), we get 


C(t, t'] v) + Fl 2 C(t,t'] v ) + / dtifj,(t — ti, v)C(£i, t 1 ; v) 

J — OO 
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where the /r(r; v) is the Fourier transform of 
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From Eq.(85), we also get the stationary solution for the dipole operator, correct to 
second order: 


aft) = 


2dQ r°° do; f d 2 k <r 3 ( 0 )ff • £^ free (k, u + k • v) 


3 i(kr— uit) 


h j-oo 2n J (2n) 2 f 2 2 [l — A (a;; v)] — cu 2 — ior^o;; v) 


where 7 ( 0 ;; v) = /x(a?; v)/2 and (P denotes the principal value) 

1,00 do/o; 2 /r(u/, v) 


A(ca;v) = -P f 
Jo 


, ( 88 ) 


(89) 


7T Q 2 U>' 2 — U 2 

are both even in u and give the second-order atomic frequency shift and decay rate. 
(They depend also on the transition dipole fj and on velocity) 

Finally, inserting Eq. ( f 88 j ) into Eq. ( [ 86 ] ) and Fourier transforming the resulting 
expression, we can write the dipole spectrum ( [81) to fourth order in the dipole 
coupling as 

Sij(u; v) = 2 h J d 2 k 0 (uj + k • v) Im R(u + k • v) oT 2kz 


k h* 

w / \ i^n^m * r \ 

x a in (w;v )—— a mj (ca;v) 


where we defined 

v ) = 


err),/if .} 2 


f2 2 (l — A(ca; v)) — u> 2 — ica 7 (a;; v) 


(90) 


(91) 


Note that this velocity-dependent polarizability differs from the 'simple damping' 
velocity-independent form given in Eq. (67), as it contains non-Markovian memory 
effects through the frequency-dependent shift A(cu; v) and damping 7 ( 0 ;; v). Using 
the symmetry in u of all involved functions, one can show that 5V, (u : v) is even in v 
and that for small velocities, it satisfies the fluctuation-dissipation relation: 

Sij(u ; v) = 2h0(u) Im 077 ( 0 ;) + 0{y 2 ) (92) 


where the imaginary part of the 'dressed' polarizability is 

Imfiy (w) = j d 2 k Im R(u) e~ 2kz a in (uj ; 0) kn ^ m affu: 0) . (93) 

Using this in Eq. ( [82] ), we obtain the quantum friction force to fourth order in the 
coupling, namely 

1 00 

F =- \ fd 2 kj e~ 2kz f do; Im R(u) 

IX J rC J 

0 

x kik* k 6{ k • v — oj) Im oy^k • v — u) . (94) 

The key observation is that the step function 9(k ■ v — u) limits the o;-integral to the 
narrow spectral range 0 < w < k • v of the anomalous Doppler effect. For small 
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velocities, we can expand the frequency-dependent functions Imf?(w) and Ini a l3 (oj) 
around uj = 0. Only the first derivatives contribute since both functions are odd in u. 
One obtains in this way 

45T),V7; 2 

F « - Ima'(0) lm#(0), (95) 

647rz' 

where d'( 0 ) is the frequency derivative of the dressed atomic polarizability given in 
Eq.(|93}, evaluated for an atom at rest (v = 0) at distance z from the surface. (For 
the full distance dependence, one has to perform the ^-integral in Eq.(|93} to obtain 
Im oiijfo) ~ 1 /z 3 .) Note that, in contrast to the prediction of the previous subsection, 
non-equilibrium fluctuation electrodynamics results in a ^''-dependence for quantum 
friction. 

A few remarks are in order. First, the next-order term proportional to v 2 in 
the expansion of the dipole power spectrum in Eq.(|92} leads to corrections to the 
frictional force proportional to v 5 . Second, we note that the result in Eq. ([95}, 
derived from the fourth-order expansion of the dipole-dipole correlation for the 
moving two-level atom, and valid in the low velocity limit, coincides with the 
result of fluctuation electrodynamics in local equilibrium Iil4l 161 [17] when the 
corresponding perturbative expression for the polarizability is used, and differs from 
the frictional power Pa in Barton's calculation by a factor of 5 [see Eq.(|24}]. Third, 
the same expression for the frictional force in Eq.([95} is obtained for the moving atom 
treated as an isotropic harmonic oscillator, a case in which exact expressions for the 
dipole-dipole correlation and a non-equilibrium fluctuation-dissipation relation are 
available Ifl8| . Finally, it is possible to show that a peculiar cancellation occurs in 
the computation of the fourth-order dipole-dipole correlator for an atom moving at 
constant velocity, which translates into an exact cancellation of terms linear in v in 
the frictional force Ii43| . 

6. Conclusion 

In summary, we have shown that the calculation of atom-surface quantum friction 
in the formulation based on perturbation theory 0 depends on how the atom is 
boosted from being initially at rest to a configuration in which it is moving at constant 
velocity parallel to the planar interface. We pointed out a subtle cancellation between 
the one-photon and part of the two-photon dissipating power. As a result the leading 
order contribution to the frictional power is quartic in velocity. Also, an alternative 
calculation of the average radiation force leads to the same conclusions, that is atom- 
surface quantum friction scales as v 3 . 
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We have reviewed recent calculations (Scheel and Buhmann [0 and Intravaia et 
al. OH) that generalize fluctuation electrodynamics for the computation of the atom- 
surface quantum friction in the non-equilibrium stationary state. They agree on the 
way the friction force is determined by the fluctuation spectrum of the dipole alone 
[Eq.(j72])], but differ in evaluating that spectrum, in particular in the low-frequency 
regime where the anomalous Doppler shift Il32l arises (ui ~ v/z). This leads in 
one case |0 to a friction force linear in v, and in the other Bl8fl to a v 3 force. To 
validate the master equation techniques behind these approaches and to resolve this 
discrepancy, it would be very interesting to extend the time-dependent perturbation 
theory pursued here and to calculate, for example, atom-field correlations in the 
stationary state. 
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Appendix A. Two-photon emission 


Appendix A.l. Leading to a force 0(v) 

The result for the frictional power in Ref.0 that turns out to scale like 0(v 2 ), arises 
from the following integral (Barton's notation P B , Eq.(4.11) of Ref.||9l with the missing 
prefactor from the erratum) 


Pb = 


ha 2 n 2 


jdmd^e- 2 ^^ 


(ki • k 2 - M’2 


(2?r) 3 

x (uq + cc 2 ) Im -R(cui) Im R( ce 2 )- 


k±k 2 

5(tt - u^)(k 2 • v) 2 


(A.l) 


(U + u; 2 ) 2 (U + c4) 2 

where, deviating from Barton's notation, the prime denotes the Doppler-shifted 
frequencies, e.g.: u\ = u — k x • v. We expand to the leading order in v and approximate 
8(Fl — cj()/(U + w' 2 ) 2 ~ <5(f2 — cui)/(f2 + w 2 ) 2 in the second line (drop the primes). 
The integrals over the wave vectors k 12 are then elementary and give 97r 2 c 2 /(16;A) 
- this frictional power is quadratic in the atomic velocity v. The (5-function fixes one 
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frequency to = Ft. The remaining frequency integral is evaluated for a narrow 
surface plasmon resonance, T <C u>s- This gives 


dee 


Tree; 


(ff + o;) 3 4u s (Fl + Co's) 3 
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<4 r 

AFIu a s 


(A.2) 


Barton gives the first term, and the second arises from the low-frequency limit 
of the surface plasmon spectral density It contributes in particular in the regime 
Ft ~ T -C ojy. Putting everything together. Barton's approach yields 
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oosiFl + c o s ) 3 


+ 


nFlujg 


(A3) 


7 


where we have marked the excited state decay rate in the short-distance limit 
[Eq.(2.11) of Ref.f9]]]. 


Appendix A.2. Leading to a force 0(v 3 ) 

Barton's result for the frictional power that turns out to scale like 0(v A ) 
the following integral (Barton's notation P A , Eq.(4.11) of Ref. 11911) 

Pa = r AAAF [ dKldK2e -^»(b k * ~ k ‘ - k ^ 2 


arises from 


167T 3 
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UiLU 2 (uJi + u 2 ) 


k 1 k 2 

5(oj'i + lu 2 ) 


|F( Wl ,w s )| 2 |F(^, Ws )| 2 (n + u»;) 2 (!J + w y 2 (A ' 4) 

where the prime denotes the non-relativistic Doppler shift: u\ = u — ki • v. Note that 
the 5-function enforces energy conservation in the frame comoving with atom: the 
pair of plasmons has zero energy there, uj( + uj 2 = 0. Since the frequencies u\, 2 in the 
laboratory frame are positive, this condition can only be satisfied if the Doppler shift 
is anomalous, for example < 0. The same condition also explains the spectrum of 
Cherenkov radiation l32l . 

To evaluate the integral ( |A.4[ >, we assume that the Doppler shift is small enough. 
More precisely, note that the exponential factor e 2(k ^+ k 3 z provides a typical range 
0(1/z) for the fc-vectors. The Cherenkov condition 0 < uq + lu 2 = (kr + k 2 ) • v then 
restricts oji j2 to the range 0(v/z) and the required approximation is |k-v| = 0(v/z ) -C 
Ft, cos- The frequency integrals then give in the leading order 

ki • k 2 ) 2 




k±k 2 


■[(kr + k 2 ) • v] 4 (A.5) 


(ki + k 2 ) • v > 0 

The restriction on the integration domain can be lifted, multiplying with f since the 
integrand is even under the transformation (k lx , k ly , k 2x: k 2y ) i->- (—k 2x , k 2y: —k lx , k iy ) 
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(v points along the x-axis). The fc-integrals then reduce to 277r 2 /(16z 10 ), and we get 
Barton's Eq.(5.4) 


P A ~ 


9 hv 4 a 2 uj 4 T 2 


512vr 


ca|z 10 


(A.6) 


Note that this expression cannot be written in terms of the (distance-dependent) 
decay rate which depends on the plasmonic mode density at the atomic resonance 14. 
The calculation above illustrates that, on the contrary, the two-plasmon emission in 
Pa is concentrated at much lower frequencies 0{v/z). In the limit Ft -C us, however, 
Eq.( A.6[ > contains exactly the same scaling compared to Eq.(|95]) of the fluctuation 
electrodynamics, and is just smaller by a factor 1/5. 

Note that it is not obvious that the radiated power P = P A + P B and the frictional 
power — F • v give the same result, as the energy taken from the atomic motion 
may also be used to excited the atom. This term, denoted by dQ/dt by Dedkov & 
Kyasov If42l . is discussed in Secs|3}|4]for different atomic trajectories. 


Appendix B. Piecewise constant acceleration 


We evaluate in this appendix the state of the atom+field system in first and second 
order of perturbation theory, using an atomic trajectory whose velocity increases 
continuously over a finite time (thin blue in Figj3]). This serves as a check of the 
general result ([38]) in the long-time limit and provides a complete list of terms that 


enter into the two-photon production rate (45). 

The one-photon amplitude c[ L \t) for the component fj, k) of the state (10) is 
proportional to 


A(e, kc a; t) — / dt j e 




(B.l) 


where r(fi) is the atomic path. This determines in the next order the amplitude of 
\ki, nf) to be proportional to 

t 

= f d£ 2 A(e,k 1 o; 1 ;£ 2 )e i (- n+a ' 2)t2 e - ik2 - r ^) . (B.2) 

— OO 

We consider a particular trajectory, namely a path with piecewise constant 
acceleration (see caption of Figj3[): 

0 t < —T 
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We define ft — k • vr, and consider the limiting case of small velocity (3 <C 1 
and 'smooth launch' |ft ± ca|r 1. For simplicity, we condense the notation into 

A(e, kca; t ) i-)- A e (t), Fl + lu ca, and set r^l. 


One-photon process. To perform the integration Eq.( |B.l[ >, we consider first the case 
that t < — 1. The integral is elementary: 


Aojt 


t < -1 : 


Ae(t) = — 


1C Jj 


(B.4) 


assuming that the coupling is switched off at the lower limit. For — 1 < t < 1, the 
ti-integral is split into —oo < 1 1 < —1, giving A e (— 1) from Eq.( |B.4[ >, and into —1 < 
ti < t, which makes a phase factor appear under the integral: e ~ ,k ' r( T) = e — i / 3 (*i+ 1 ) 2 / 4 . 
Since this phase is < (3 <C 1, we expand this exponential and find 
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When the acceleration is finished, we thus get 
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(B.6) 


Finally, for larger times, the branch of the path with a constant velocity contributes 
between 1 < t x < t, again an elementary integral: 

g i {u-p)t _ e i(w-/ 3 ) 

t > 1 : Ae(t) = A e ( 1) + 


i ca 


P) 


* 4 -e(l) + T. 




piua piua piua 

— + (3— + if3—. 


m ■ ■ - - - 2 ■ (B.7) 

i(ca — p) ica ca ca 2 

We have expanded to first order in [3 all terms except the one where [3t appears in 
the exponent, because we shall be interested in the long-time limit. Note the three 
cancellations with A e {l) so that we get 


t > 1 : A e {t) 


-p)t 


e‘ v ~ (3(e u 


i(ca — /3) 2ca 3 i(ca — (3) 

Putting everything together and restoring the physical units, we get 

gi(f 2 +o;i)t 

i(f2 + cai) 


Q i{u-p)t i/3 silica m 

+ 4-• (B.8) 


ca 2 ca 


*4 e (kicai; t) = < 


e i(Q+o iki . v sin (Q + Wl ) r 
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. i(ff + cai —/3i) (fl + cai) 2 (f2 + cai)r 


for t < —t 

(B.9) 

for t > r 
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Note that the second line is exactly of the form put forward on general grounds in 
Eq.([38]) of the main text. The first term corresponds to the 'adiabatic limit' where the 
atomic velocity is taken at its final value. It is independent of the duration r of the 
acceleration. In the second term, the sine function (last fraction) reduces to unity for a 
sudden acceleration (limit r —> 0). Any finite value of r decreases this amplitude, and 
effectively suppresses it when the atom is smoothly accelerated, i.e., (fl + ui)t 1. 


Two-photon process. Its amplitude is given by integrating the one-photon amplitude 
A(e, kicui; t 2 ) once again [Eq.( B.2[ >]. We use the previous expression (B.9l and Eq.(B.7). 
Two small parameters appear /?i :2 = k 12 • vr that we consider small and of the same 
order. 

We begin for t < — 1 with an elementary integral 

e i(pl+W2 )t 

(n + Unicoi + U 12 ) 

For — 1 < t < 1, the first-order expansion in /3i, /3 2 yields: 

- 1 < t < 1 : M(t) « M(- 1) 


(B.10) 


-1 


pi(ui+W 2 )t 2 / ^ 

dt 2 — - -fl-^(t 2 + l) 2 

i(ff + cai) V 4 v ’ , 


/A [ d U e i(wi+W2)t2 { ^ + ; + ^ 1 

J \A(n + uj 


-1 

t 


+ Pi Jdt 2 e i ( _f2 +^2 ) t2 

-1 


\4(n + o;i) ^(ff + cai) 2 + 2(ff + cai) 3 J 

—i(Q+cdi) 

(B.ll) 


2(11 +a;i ) 3 ’ 


where the last two lines arise from Eq.( |B.7| . The second line is an integral analogous 
to Eq.( |B.5[ >, and the result partially cancels with M{— 1). The other integrals are just 
a bit tedious to work out and eventually yield the cumbersome expression 

e i(u)l+a) 2 )i 


M(t) 


(12 + CUi)(cUi + UJ 2 ) 

Pi + Pi c i (u 1+ u 2 )t j (t + !) 

i(f2 + cai) 


. t + 1 1 

\ 4(cei + oo 2 ) 1 2(w 1 + ca 2 ) 2 + 2(cai + ca 2 ) 3 


Pi + P 2 e-^ 1 ^ 

i(II + tai) 2(c+ + ca 2 ) 3 

i/3i - i(c+ + u 2 )(t + 1)) - 


2(n + w 1 ) 2 
Pi 


(cai + cu 2 ) 2 


g i(o;i+W2 )t _ g-i+l+^2) 


2(Fl + tai) 3 i(cui + c of) 
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(B.12) 


2 (ft + cai ) 3 i(—ft + o^) 

We finally get to the physically interesting case of late times where Eq.( |B.8[ ) can be 
used and the integrals become again elementary 


L 

t > 1 : M(t ) = M( 1) + J 


d U 


Qi(u 1 - 0 x+ui 2 - 0 2 )t 2 

i(fl + UJ\ — /3f) 


r dt cl[ -n+u 2 - 02 )t 2 ^i sin(tt + a>i) 
(ft + cai ) 3 


= M{ 1)- 


e i(u} 1 -0i+uj 2 -0 2 )t _ q1(u)i—0i+lc 2 —0 2 ) 


+ 


(ft + CUi - Pi)(ui - P1+UJ2- P2) 
e i(-o+^2 -0 2 )t _ e i(-n+u 2 -0 2 ) ^ sin (-Q + 


(B.13) 


—ft + UJ2 — P 2 (ft + cei ) 3 

This can be written as a sum M A {t)-\- M Me whose terms we discuss separately 

now. 

The 'adiabatic amplitude' appears in the first line of Eq.( |B.13[ > 

e iK+< 4 )t ( 2 ft + u[ + u' 2 ) eM+wji 


M A (t) = 


->■ 


(B.14) 


(ft T cu()(cu( T cu(>) (ft T cu()(fl T T ce^) 

where the notation u[ = — was used. The term after the arrow (—>) in this formula 

gives the total amplitude after symmetrizing the quantum numbers k^i and k 2 o ;2 of 
the two plasmons. This expression is independent of r because M must be a squared 
time by definition. It is identical to the term featuring e lAi in Eq.(4.6) of Ref-H, the 
one that leads for t — > 00 to the 5-function 8{u[ + ujf) and the power P A . 

The second line of Eq.( |B.13 > contains the other time-dependent term: 

( ki . v ) sin (Q + Wl ) r 


Ms{t) — 


(B.15) 


-Fl + u '2 (ft + ca!) 3 r 

where we have restored r. To this order in v, the limit r —>• 0 recovers the term 
proportional to eK-^+Ao)* in Eq.(4.6) of Ref.|9j]. We recall that this term leads to the 
5-function 5{u' 2 — ft) and the power P B scaling with v 2 , together with its exchange- 
symmetric partner. Note that this term, up to the first factor, is exactly given by the 
second (constant) piece of the one-plasmon amplitude A e ( ki^p t) in Eq.( B.9[ >. Hence 
the sine function reduces its contribution if (ft + ujf)T 1. 

The remaining term Me collects all terms independent of t in Eq.( |B.13| >. Their 
expansion for small /?i i2 is tedious and leads to 


M c = 


Pi 


+ 


P2 \ sin(a;i + u 2 ) 


(bJ\ + CO2 ) 2 ( fl — SJ2 ft + CCi J (bJ\ + CU 2 ) 


+ 0{p 2 ) . (B.16) 
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We add the corresponding expression under plasmon exchange (1 2 ) and get in 

physical units 


M 


2ft 


C,sym 


ki ■v k ? •v 

+ 


sin(cai + lu 2 )t 


(B.l 7) 


(oJl + UJ 2 ) 2 { H 2 — (jj\ ' fl 2 — Ulj } (cui + (jJ 2 )t 

It is straightforward to check that this is equal to the small-v expansion of the 
constant terms in M(t) as given in Eq.(4.6) of Ref. |9| . 

To summarize, we have extended the calculation of the complete two-photon 
amplitude to an atomic path with an acceleration phase of duration 2 r. The two- 
plasmon power called Pa, scaling with 0(v 4 ) does not depend on r and is unchanged, 
at least for small velocities. The power P Br scaling with 0(v 2 ), depends on r and 
becomes suppressed when the launch duration r is larger than the atomic period 
1/0. We have confirmed the argument given earlier that this 0(v 2 ) power can 
be computed from the first-order transition amplitude: it is proportional to the 
probability of exciting the atom in a non-adiabatic way This means: the acceleration 
has lead to an amplitude shift in the 'Lamb cloud' of virtual photons surrounding 
the atom. We may say that these photons have become 'real' because their amplitude 
differs from the adiabatic value. 

Technical note: To get a probability amplitude that increases linearly with t, we need 


-Avt 


lim 

fc—> OO 


V 


= 2tx\5(u) , 


lim 

£—>•00 


-Aut 


V 


= 2wt5(v) . 


If the '— V is rather a complex function a(v), we may evaluate 


e ivt - 1 + 1 - a(v) 


v 


—» 2Trt5(v) + 47 t 5 ( z /) 


Im a{v) 


+ 


|1 — a(u) 


v 


(B.18) 


(B.19) 


where the last two terms do not increase with t (if the final integral converges, of 
course). Hence they drop out when a transition rate is calculated. 


Appendix C. Evaluation of the average force 

Appendix C.l. Second order: exponentially small 


The second-order term of the average force (271 is given by 

F t2 /f) = 2Re^fJd 3 n ( g , vac|F(f)|rf, «:)ci 1) (f) 
v 

= 2 V Re ^/ d3Kk|,? ' £|2W2 (ii + J--iA) ’ 


(C.l) 

(C.2) 
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where we have used the matrix element (111 and the amplitude c^(t) [Eq.(14t], We 


consider in this appendix only the long-time limit where r(£) = vt. Summing over 
the excited states and taking the real part, one gets: 

O /~/2 p 

F (2) (t) = -/ d 3 K k k e~ 2kz Im R(u)6{n + J) (C.3) 

7r J 

which is nothing but Eq.([28]>. The resonance condition Fl+u' = 0 can only be satisfied 
for large k = 0(Fl/v), making this contribution exponentially small in v. 


Appendix C.2. Fourth order, via vacuum 

We continue with the fourth-order part involving the mixed amplitude c^*(t)c^\t) 
in Eq.(27). This product can be combined with the last line of Eq.([2l| for (t) where 
we recognize the expression for c^(f). The sum yields a force (subscript 0 for 'going 
via zero-photon sector') 

Fo 4) (t) = 2R e Il/d 3 Ac(^,vac|F(f)|7f,/c)cS 1) (f) 

* {-*- *£*-'* }■ <C4) 
Here, we recognize one term, — 7s f F (2 ) (t ), quadratic in exponentially small parts, that 
translates the loss of probability in the ground state. For the other piece, we use the 
identity 


Re k 


= Re { (V2 + 


~ 2' VvRe fi + w'-iA 


SE ri 


V v Im 


h + a/ — iA 

to identify ground-state decay rate and level shift from Eq.([T7|) 

Fo 4) (t) = - 2Re^/d 3 ^ {g, vac|F(t) \ff, k)c \^ (t) 


(C.5) 


= 7 j d 3 K\rj • £| 2 |0 K | 2 Re 


1 


£1 + uj f — iA 


V 

= -V v ( 7 g 6E g ). (C.6) 

Again, this is an exponentially small term. For its interpretation, one may think about 
the adiabatically stored energy in the Lamb-shifted ground state 5E g . 
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r-r-rl / OX 

The final piece for the force arises from that part of c\ it) that goes via the two- 
photon sector [first and second lines of Eq.pT])]. We have to add the mixed term 
from the coherence between the one- and two-photon sectors. The two contributions 
are denoted F^ 03 ' and F^ 12 ^ and are handled separately (subscript 2 for 'going via 
two-photon sector'): 

i d 3 f* K e i(n+w '^ 


p (4) [03] 


= 2R ejfj d 3 K(^,vac|F(t)|/7,K)|4 
x J d 3 Ki 


3 (14 + oj' - iA) 

(v ■ £i)|0 Ki | 2 (£i • k)* ( 1 


-\- uj f — iA) 


= Re E /d 3 «:d 3 Ki \ki ■ k\ 2 \f K \ 2 \f 


+ cjq — iA + uJ — iA y 

2 _ik_ 

h ' 11 (14 + u' — iA) {u'i + u' — iA) 


x 


{; 


1 


1 - } . 

/-at 


(C.7) 


■ 14 T — iA 14 T oj 1 — iA 

To proceed, we neglect exponentially small terms arising from the 5-functions 5(14 + 
u') and drop the —iA in the corresponding non-resonant denominators. The only 
term that remains is (re-labeling and Bose symmetrizing the integrand) 


^(4) [03] _ 
2 — 


7rd 4 

If 


X 


?/ 

V 

ki 


d 3 Kid 3 K 2 |A:i • fc 2 | 2 |0 Kl | 2 |0 K 2 | 2 <5(a ;' 1 + a/,) 
k2 'i 2f2 -f- C-f- Co^2 


7 + 




• c -P ^2 ^ ^i)-j - 

Finally, we have to add the triple integral (symmetry factor 1/2 cancels) 

f ( 4 )[1 A = Rejfl d 3 Kd 3 Ki d 3 K 2 (rf, n\F(t)\g, {t)c [ 2 \t) 


(C. 8 ) 


(C.9) 


Insert the matrix element (12) with its two Bose-symmetric terms, the amplitudes ( [14} 
, exploit the 5-functions 5(k — ki j2 ), to get 


p,(4)[12] 


= ^3 Re / d 3 Kid 3 K 2 |fci • k 2 \ 2 \<j) K1 \ 2 \f 
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(_i_(. 

1 -)- cjo — iA ' 1 
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(C.10) 

'12 + wi-iA n + tu'z-iXn v ' 

Neglecting again the imaginary part in the non-resonant denominators we observe 
that the same structure as Eq.( |C. 8 } emerges, up to a switch k, <h- k 2 in the photon 
momenta. By symmetry, we expect that F and v are parallel and find for the 





















Friction forces on atoms after acceleration 
projection 

v . (Fi 4)[03] + F^ 4)[12] ) ~ - 


41 


7 rd 4 

~w 
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2 Fl 4" tufj \ 2 


x v ■ (k x + k 


(C.ll) 


- (n + cuj)(f2 + u'f) ' 

Under the 5(o;j + t o'f), we may replace v ■ (kx + k 2 ) ha uj 1 + u 2 , and we recover the 
structure of the two-photon emission 

nd j d3 «i d3fi: 2 \ki • fc 2 | 2 |0 Kl | J |0 K2 | 2 5(cax +w 2 ) 
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4 U 2 


v • (F 2 4 ) [03] + F 2 4)[12 ^) 
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